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Abstract 

We study some rare decays of Bg meson governed by the quark level transitions 6 — )• s, in the 
fourth generation model popularly known as SM4. Recently it has been shown that SM4, which is a 
simple extension of the SMS, can successfully explain several anomalies observed in the CP violation 
parameters of B and Bg mesons. We find that in this model due to the additional contributions 
coming from the heavy t' quark in the loop, the branching ratios and other observables in rare Bg 
decays deviate significantly from their SM values. Some of these modes are within the reach of 
LHCb experiment and search for such channels are strongly argued. 

PACS numbers: 13.20.He, 13.25.Hw, 12.60.-1, ll.30.Er 
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I. INTRODUCTION 



The spectacular performance of the two asymmetric B factories Belle and Babar provided 
us an unique opportunity to understand the origin of CP violation in a very precise way. 
Although, the results from the B factories do not provide us any clear evidence of new 
physics, but there are few cases observed in the last few years, which have 2-3 a deviations 
from their corresponding SM expectations [1]. For example, the difference between the direct 
CP asymmetry parameters between B~ 7r^K~ and B^ — t- 7c~^K~, which is expected to 
be negligibly small in the SM, but found to be nearly 15%. The measurement of mixing- 
induced CP asymmetry in several 6 — )■ s penguin decays is not found to be same as that of 
Bd — J/ipKg. Recently, a very largish CP asymmetry has been observed by the CDF and 
DO collaborations [2], Sj in the tagged analysis of Bg — t- J/ip(j) with value S^^ G [0.24, 1.36]. 
Within the SM this asymmetry is expected to be vanishingly small, which basically comes 
from Bg — Bg mixing phase. It should be noted that all these deviations are associated with 
the flavour changing neutral current (FCNC) transitions 6 — )■ s. It is well known that the 
FCNC decays are forbidden at the tree level in the standard model (SM) and therefore play 
a very crucial role to look for the possible existence of new physics (NP). 

In this paper we would like to study some rare decays of Bg meson involving 6 — )■ s tran- 
sitions. The study of Bg meson has attracted significant attention in recent times because 
huge number of Bg mesons are expected to be produced in the currently running LHCb ex- 
periment, which opened up the possibility to study Bg meson with high statistical precision. 
These studies will not only play a dominant role to corroborate the results of B mesons 
but also look for possible hints of new physics. Here we will consider the decay channels 
Bg —7- 07r, Bg —7- 07, Bg —7- 77 and Bg — )■ fi^fi^'y which are highly suppressed in the SM. We 
intend to analyze these decay channels both in the SM and in the fourth quark generation 
model ^, usually known as SM4. SM4 is a simple extension of the standard model with 
three generations (SM3) with the additional up-type (t') and down-type {d') quarks. It has 
been shown in Ref. {s], that the addition of a fourth family of quarks with mf in the range 
(400-600) GeV provides a simple explanation for the several deviations, that have been ob- 
served involving CP asymmetries in the B, Bg decays. The implications of fourth generation 
in various B decays are discussed in j6-10|. The experimental search for fourth generation 
quarks has also received significant attention recently due to the operation of Large Hadron 
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Collider. The CMS collaboration put a lower bound on the mass of t' as m^/ > 450 GeV 



and exclude the 6'-quark mass in the region 255 GeV < my < 361 GeV at 95% C.L. 12|. 



The paper is organized as follows. In section II we discuss the non-leptonic decay process 
Bg —7- (j)7i. The radiative decays Bg — t- 07 and Bg — 77 are discussed in Sections III and IV. 
The process Bg — )• fi^^~'y is presented in Section V and Section VI contains the Conclusion. 



II. Bg 



bn PROCESS 



In this section we will discuss the non-leptonic decay mode Bg — t- (j)TT which receives 
dominant contribution from electroweak penguins b — > sqq (g = u,d), as the QCD penguins 
are OZI suppressed, and the color-suppressed tree contribution b — )■ uus is doubly Cabibbo 
suppressed. Therefore, this process is expected to be highly suppressed in the SM and hence 
serves as a suitable place to search for new physics. This decay mode has been studied in the 



131 ] and in the model with non-universal Z' boson 



SM using QCD factorization approach 

The relevant effective Hamiltonian describing this process is given by 



SM 



V2 
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(1) 



1=1,2 i=3 

where Cj(/i)'s are the Wilson coefficients evaluated at the 6-quark mass scale, Oi,2 are the tree 
level current-current operators, Os-e are the QCD and O^^iq are the electroweak penguin 
operators. 

Here we will use the QCD factorization approach to evaluate the hadronic matrix ele- 
ments as discussed in |l^. The matrix elements describing the Bg (p transition can be 



parameterized in terms of various form factors [17[ as 



(</)(p',e)|s7^(l - -f^)b\Bg{p)) = -iel{mB,+m^)Ai{q^) 

+ ^{V + V),{e* ■ q) ^^^^ + iq,{e* ■ q)'^-^{A,{q') - Ao{q')) 



*u p la 



2V{q') 



(2) 



where p and p' are the momenta of Bg and cf) mesons, q = p — p' is the momentum transfer, 
Ai_^{q^) and V{q^) are various form factors describing the transition process. Using the 
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decay constant of 7r° meson as 



71 " ^7! ' 

one can obtain the transition amplitude for the process 



(3) 



(4) 



where Ag = VqbV*^. The parameters a^'s are related to the Wilson coefficients Cj's and the 
corresponding expressions can be found in Ref. 16 1. 
The corresponding decay width is given as 



e* ■ g 



(5) 



where Pcm is the center of mass momentum of the outgoing particles. 

Now we discuss about the CP violating observables for this process. To obtain these 
observables, we can symbolically represent the amplitude (jl]) as 



A{B, ^ 07r) = (e* ■ q)[\uA^ - XtAt] = -XtA^e* ■ q) 



1 — r a e 



(6) 



where a = |A„/At|, —7 is the weak phase of Vub, (7r + /3s) is the weak phase of Xt, r = \Au/At\, 
and 6 is the relative strong phases between At and A^. From the above amplitude, the direct 
and mixing induced CP asymmetry parameters can be obtained as 

2ra sin 6 sm{(3s + 7) 



A^ 



1 + {ray + 2ra cos S cos{f3s + 7) 
2ra cos 6 sin(/3s + 7) + {ra)"^ sin(2/3s + 27) 
1 + {ray + 2ra cos 6 cos{/3s + 7) 



(7) 



For numerical evaluation, we use the particles masses, lifetime of Bg meson from 18|. For 
the CKM elements we use the Wolfenstein parametrization with the values of the parameters 
as A = 0.2253 ± 0.0007, A = 0.808l|^:[!?^, p = 0.132l°:°?2, fj = 0.341 ± 0.013. The parameters 
of QCD factorization approach and the value of the form factor used An"^'^ = 0.32 ± 0.01 

n 

are taken from [13|. 

With these inputs we obtain the branching ratio for this process as 



Bt{Bs (f)7i) = (1.26 ± 0.32) X 10 
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(8) 



which is consistent with the prediction of 13|, [14 1 . 



The CP violating observables are found to be 



-0.23 



A, 



0.1 



(9) 



Our predicted direct CP asymmetry is lower than the prediction of 13|. This difference arises 
mainly because the sub-leading power corrections to the color suppressed tree amplitude 02 
has been included in Ref. 13|, which introduces a large strong phase. 



Now we will analyze this process in the fourth generation model. In the presence of a 
sequential fourth generation there will be additional contributions due to the t' quark in the 
loop diagrams. Furthermore, due to the additional fourth generation there will be mixing 
between the b' quark the three down-type quarks of the standard model and the resulting 
mixing matrix will become a 4 x 4 matrix (Vckma) and the unitarity condition becomes 
+ K + K + K' =0, where Xg = VghV*^. The parametrization of this unitary matrix 
requires six mixing angles and three phases. The existence of the two extra phases provides 
the possibility of extra source of CP violation 19|. In the presence of fourth generation 
there will be additional contribution both to the Bg — t- (pir decay amplitude as well as to 
the Bg — Bg mixing phenomenon. However, since the new physics contribution to B^ — Bg 
mixing amplitude due to fourth generation model has been discussed in Ref. |9!], we will 
simply quote the results from there. 

Now we will consider the additional contribution to the decay amplitude due to the fourth 
quark generation model. In this model the new contributions are due to the t' quark in the 
penguin loops. Thus, the modified Hamiltonian becomes 

10 10 



(10) 



XuiCiO, + C2O2) -XtJ2 c^o, - Xt> J2 C[Oi 

where Cfs are the effective Wilson coefficients due to t' quark in the loop. To find the 
new contribution due to the fourth generation effect, first we have to evaluate the new 
Wilson coefficients C[. The values of these coefficients at the Mw scale can be obtained 
from the corresponding contribution from t quark by replacing the mass of t quark by t' 

lese values can then be evolved to the mi, scale 
15| . Thus, the obtained values of C[^j_iQ{mh) for 



mass in the Inami Lim functions [20]. T 
using the renormalization group equation 



two representative set of values i.e., mj/ = 400 and 500 GeV are as presented in Table-I. 
Thus, in the presence of fourth generation model, one can obtain the transition amplitude 
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TABLE I: Numerical values of the Wilson coefficients C'^ for mt' = 400 and 500 GeV. 



t' mass 




C's 




r" 

'-'10 


mf/=400 GeV 
mt, =500 GeV 


4.453 X 10^3 
7.311 X 10-3 


2.115 X 10^3 
3.199 X 10-3 


- 0.029 

- 0.041 


0.006 
0.009 



for Bs —7- </)7r process from Eq. ( JTOl) . which can be symbohcally represented as 



(11) 



where b = \Xt'/Xt\, r' = \At//At\, and 6i is the relative strong phases between Afi and At. 
From the above amphtude, the CP averaged branching ratio, direct and mixing induced CP 
asymmetry parameters can be obtained as 



Br = Br^^^ X 



Ad,n 



Y 

X ' 



Sdm — 



X' 



(12) 



with 



X = 1 + {raY + {r'by + 2ra cos 6 cos(/3s + 7) + 2r'b cos 61 cos(0s — f3s) 

+ 2rr'abcos{(f)s + 7) cos(5 + 61) , 

Y = 2ra sin 6 sin(/3s + 7) + 2r'b sin 61 sin(0s — f3s) + 2rr'ab sin(05 + 7) sin((5 + 5i) , 

Z = sin 29 + 2ra cos 5 sin(/3s + 7 + 29) - 2r'b cos 5i sin(0^ - /3s - 26^) 

+ r^a^ sin(2/3, + 27 + 2^) - r'^fe^ sin(20, - 2/3, - 20) 

- 2rr'a6cos((5 + 5i)sin(0, -2/3, -7-20). (13) 



In Eq. ( ]T3ll . 26* is the additional contribution to the Bg — Bs mixing phase in the fourth 
generation and the expression for it can be found in Ref. 9|. 

For numerical evaluation using the values of the new Wilson coefficients as presented 
in Table-I, we obtain r ^ 7.79, 5 ^ 25.9°, r' = 3.48 (5.03), and 61 ^ -0.1° (-0.1°) 
for rrit' = 400 (500) GeV. For the new CKM elements Xt', we use the allowed range of 
I I = (0.08 - 1.4) X 10-2 [(0.06 - 0.9) x lO'^] and 0, = (0 ^ 80)° [0, = (0 ^ 80)°] 
for nit' = 400 GeV [500 GeV], extracted using the available observables mediated through 
6 — )■ s transitions . Now varying and 0, in their allowed ranges, we show the variation 



of branching ratio in the left panel of Figure-1 and the correlation plot between the CP 
violating parameters in the right panel. From the figure it can be seen that the branching 
ratio is significantly enhanced from its SM value and large mixing- induced CP violation 
(S^n) could be possible for this decay mode in the fourth generation model. However, the 
direct CP asymmetry does not deviate much from the corresponding SM value. It should 
also be noted that the branching ratio decreases slowly with the increase of t' mass. However, 
there is no significant rrit' dependence of the CP violating observables. 




FIG. 1: Variation of the CP averaged Branching ratio in units of 10^ (left panel) and the Correla- 
tion plot between the mixing induced CP asymmetry 8^-,^ and the direct CP asymmetry parameter 
Atf)Tr (right panel) for the Bg — )• (pir process. The red (blue) regions correspond to mt'=iOO GeV 
(500 GeV) 



III. B, ^ (Pj 

Here we will consider the decay channel Bg — )■ </)7 which is induced by the quark level 
transition b — )■ s'j. This mode is the strange counterpart of the B — )■ K*'~f, which is very 
clean to analyze. Compared to Bd meson the new elements of Bg mesons are the small 
mixing phase (pg and the large width difference AF^ of the Bg meson. The branching ratio 



of this mode is recently reported by the Belle collaboration 21 1 



Br(5, ^ 07) = (57lJ^(stat)l}2(syst)) x 10^^ . (14) 
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In the standard model the CP averaged branching ratio of this mode is predicted to be 22 1 

Bi{Bs 07) = (39.4 ± 10.7 ± 5.3) x 10"^ . (15) 



Although, the SM prediction seems to be consistent with the observed value, but the presence 
of large experimental uncertainties makes it difficult to infer/rule out the presence of new 
physics from this mode. 

The transition process 6 — )■ 57 can be described by the dipole type effective Hamiltonian 



which is given as [23 1 



4C 

n,ff = --^XtC-r{mk)Oj, (16) 
where C7 is the Wilson coefficient and O7 is the electromagnetic dipole operator given as 

Oj = ^^F^,[m,m^''{l + 75)6 + mssa'^^il - 75)6] (17) 



The expression for calculating the Wilson coefficient Cj{^) is given in 23|. 



The matrix elements of the various hadronic currents between initial and the final 



meson, which are parameterized in terms of various form factors as 17| 



(0(p',e|sor^,g^(l + 75)&|5,) = «e^,,.e*Vp"^2Ti(g2) 

+ Uq^WMl^ - ml) - {e* ■ q){p + p'),] , (18) 

with Ti(0) = T2(0) and q = p — p'. With these definition of form factors, one can obtain the 
corresponding decay width as 

^ <^7) = ^|V^t.^.:nC^r^^Pm>ljTi(0)p (1 - • (19) 



Using the value of the form factor Ti(0) = 0.349 ± 0.033 (l7|, ^7(772^) = -0.31, and the 



values of the other parameters as discussed in section II, we obtain the branching ratio as 

Br(5, ^ 07) = (39.9 ± 12.3) x 10"^ . (20) 

As is well known the rare radiative decays of B mesons are particularly sensitive to the 
contributions from new physics. The V — A structure of the weak interactions can be tested 
in FCNC decays of the type b — > {s, d)'~f, since the emitted photon is predominantly left 
handed. The crucial point is that the leading operator O7 ~ scr^^F^'^b^n) necessitates a 
helicity flip on the external quark legs, which introduces a natural hierarchy between the 
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left and right handed component of the order irid^s/fnb- However it is difficult to measure 
the helicity of photon directly. It was pointed out long back that the time dependent CP 
asymmetry is an indirect measure of the photon helicity 2J], since it is caused by the 
interference of left and right handed helicity amplitudes. The final state in Bg — t- 07 is 
not a pure CP eigenstates. Rather in the SM they consist of equal mixture of positive and 
negative eigenvalues. Thus, due to an almost complete cancelation between positive and 
negative CP eigenstates, the asymmetries in 6 —t- 37 is very small. They are given by ms/mb 
where the quark masses are current quark masses. 



The normalized CP asymmetry for the Bg — )■ (pj is defined as follows 



25| 



r(B. ^ ^,7) + r(s, ^ ,#,7) ' ^ ' 

where the left and right handed photon contributions are added incoherently i.e., T{Bs 
07) = T{Bs — 07l) + T{Bs —7- 07_r). It is well known that, the neutral mesons exhibit 
the time dependent CP asymmetry through mixing, i.e., if the particle and the antiparticle 
decay into a common final state /. In 5^ — 07 this accounts to 

Bs ^ 07L(ij) Bs (22) 
With \q/p\ = 1, the CP asymmetry assumes the following generic time dependent form 

Acpit) = Ar y ^ — (23) 

^' cosh^ -iysinhAr,t2 ^ ^ 

In terms of the left and right handed amplitudes 

^L(R) = ^(5, ^ 07L(ij)), Al{R) = A{Bs ^ (j)-fL{R)), (24) 

the form of the observables C, S and H can be found as 

C 







Ar\ 




{\Al 


P + \Ar 










P + 


Al\ 


' + \An\ 


2 



2lm[l{ALAl + ArA*j,)] 

I /I 19 I I ^ 19 I I 19 I I 4 



Al\' + \Ar\^ + \Al\^ + \Ar^^ 



2Re[l{ALAl + ARA*j,)] 
\Al\' + \Ar\^ + \Al\' + \Ar\^' ^ ^ 



In the standard model the leading operator O7, which allows the Bs{Bs) meson to decay 
predominantly into a left (right) handed photon whereas Bs{Bs) meson decaying into the 
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left (right) handed photon suppressed by an ms/rrih chirahty factor. Due to the interference 
between mixing and decay in Bs — )■ 07, a single weak decay amplitude proportional to At is 
exactly canceled by the mixing phase and hence one can obtain S^jyy = and if^^ = 2'ms/mb 



25|. 



The situation can be significantly modified in certain models beyond the standard model 
by new terms in the decay amplitudes and also by the new contribution to the Bg — Bg 
mixing. In this section we will study the effect of fourth quark generation on the various 
decay observables. In the presence of fourth generation, the Wilson coefficients Cj will be 
modified due to the new contributions arising from the virtual t' quark in the loop. Thus, 
these modified coefficients can be represented as 

Cr(/x) = CM + ^C;(/.). (26) 

The new coefficients Cj can be calculated at the M^y scale by replacing the t-quark mass 
by rrit' in the loop functions. These coefficients then to be evolved to the b scale using 
the renormalization group equation as discussed in The values of the new Wilson 

coefficients at the scale for m^/ = 400 GeV is given by Cy(mf,) = —0.375. 

Thus, including the new physics contribution due to fourth generation effect the branching 
ratio can be obtained from Eq. f ll9p by replacing C7 by C^"* and the CP violating parameters 
are given as 

/ -C| sin 29 + 2aC7C; sin(0, - - 29) + a^C',^ sin 2(0, - (3^ - 9) \ 
~ mfe V C7 + r^C? + '^aC-jC'^ cos(0, - /?,) j ^ K') 

_ m, fC^cos29 + a^C'^^ cos 2{(j)g - /3, - 0) + 2aC7C^cos(0, - /3, - 2^)\ 

where 29 is the new contribution to B^ — Bg mixing phase due to fourth generation. Now 
varying Xf between (0.08 — 1.4) x 10^^ and 0^ between (0 — 80)° we show in Figure-2, the 
CP averaged branching ratio (left panel) and the correlation plot between the CP violating 
observables (right panel). From the figure it can be seen that small but nonzero CP violating 
observables could be possible in the fourth generation model, while the branching ratio still 
consistent with the observed value. Furthermore, in this case also the branching ratio 
decreases with the increase of t'-mass. 
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0.002 0.004 0.006 0.008 0.01 0.012 0.014 -0.01 -0.005 0.005 0.01 0.015 0.02 0.025 0.03 

FIG. 2: Variation of the CP averaged Branching ratio (left panel) and the Correlation plot between 
the CP violating observables S^^/ and H^^^ (right panel) for the Bg — )• (f)-y process. The horizontal 
blue line on the left panel is the central value of the measured branching ratio whereas the green 
lines represent the corresponding 1-sigma range. 

IV. B, ^ 77 

Now we will discuss the decay process Bg — )■ 77. At the quark level this process is similar 
to 6 — 7- 57. Up to the correction of order the effective Hamiltonian for b — ?■ 577 

at scales /ib = 0{mb) is identical to the one for b — > 57. The — )• 77 process has been 



studied extensively in the SM and in various new physics scenarios 
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experimental limit on the decay -Bs — ?■ 77 is 21 1 



43 1|. The present 



Br{Bg -> 77) < 8.7 X 10^*^ (90% C.L.). (29) 

We expect with the continuous accumulation of the experimental data, the situation will 
improve and the branching ratio will be more precise. 

The effective Hamiltonian for this process is given by Eq. ( JT6l) . To calculate the decay 



amplitude for this process one may follow the procedure discussed in Ref. [29|. In order 
to calculate the matrix element of Eq. f lT6|) for Bgips) ^ 7(^1)7(^2)? one can work in the 
weak binding approximation and assume that both the b and s quarks are at rest in the Bs 
meson and the b quark carries most of the meson energy and its four velocity can be treated 
as equal to that of Bg. Hence one may write b quark momentum as pb = iTibV, where v is 
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the common four velocity of b and Bg. Thus we have 

Pb- ki = rubV ■ ki = ^rribmB, = Pb ■ h 

Ps-h^{p-h- h) ■ h = -^^(mfi^ - nib) = • ^2 (30) 
The amphtude for Bg 77 can be computed using the following matrix elements 

{0\s^^l5bmpB)) = ifB.PB^ 

{0\s'y,b\Bg) = iJbMb. (31) 

Thus, one can obtain the total amplitude for this process containing CP even and CP 
odd parts as 

A{Bs ^ 77) = M+F^.Ff^" + iM-F^^F'^'' , (32) 

with 

M-^ = -^-^^MZf^„,,y^,V* ^BmbKiml) + ^) , (33) 

and 

M- = ^^^^fB.rriBytbV* mBA,J{ml) + m,SL(m^) + ^ j , (34) 

where A = ms^ —mb- The parameters AqS are related to the Wilson coefficients Q's, which 
are evaluated at the rub scale as 

Au = {C,-C^)N,+ {C,-C,) 

A, = (Ci + C3-C5)iVc + C2 + C4-C6 

As = =Ab = -^[{C^ + C^-C,)N,+ {C^ + C^-C^)] 
B = C^~{C^N, + C^) (35) 
The functions J{m^), K{m^) and L{w?) are defined as 

Kim") = 4/n(m2)-/oo(m2) 

L{m^) = 7oo(m') , (36) 

with 
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Thus, one can obtain the decay width of Bg — )■ 77 is 



T{Bs ^ 77) 



IGvr 



(|M+|2 + |M"|2) 



(38) 



To obtain the numerical results we use the parameters as presented in section II. Thus, 
we obtain the branching ratio as 



Bt{B, -> 77) = (1.8 ± 0.4) X 10"^ 



(39) 



2l|. 



which is lower than the present experimental upper bound 

In the sequential fourth generation model there exist additional contribution to 6 — )■ S7 
induced by the 4th generation up type quarks t'. The new Wilson coefficients can be obtained 
from those of their t counter parts by replacing the mass of t quark by t' at the Mw scale, 
which is then evolved to the mi, scale using the renormalization group approach. As discussed 
in the previous section the values of the new Wilson coefficients at the rrib scale for m^/ = 400 
GeV is given by C'-j{mij) = —0.375. At the scale rrih, the modified Wilson coefficient of the 
dipole operator becomes 



Cr(m,) =C7(m,) + 



Vt'bVtr 



(40) 



Now varying Xf between (0.08 — 1.4) x 10"^ and (pg between (0 — 80)° we show in Figure-3, 
the branching ratio for Bs — t- 77 process. From the figure it can be the branching ratio can 
be enhanced from its SM value, but the enhancement is not so significant. 



1^ 
o 



X 
en 




0.006 0.008 



0.012 0.014 



FIG. 3: Variation of the Branching ratio for the process Bg — )• 77 process. 
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V. 5° ^ /i+^"7 PROCESS 



Now let us consider the radiative di-leptonic decay modes Bg — )■ /i^/i~7, which are also 
very sensitive to the existence of new physics beyond the SM. Due to the presence of the 
photon in the final state, this decay mode is free from helicity suppression, but it is further 
suppressed by a factor of a with respect to the pure leptonic Bg — ^ fi^fi^ process. However, 
in spite of this a suppression, the radiative leptonic decay B^ — t- fi^fi^j, has comparable 
decay rate as that of purely leptonic ones. 



The effective Hamiltonian describing this process Bg fi~^fi is 15 1 



n 



eff 



Gpa 

V27T 



VtbV* 



CI" {s 7^ Pl h){l Y + Cio {s 7^ Pl b){l Y 75 



2C7 nif, 




sia^,q''PRb){lY I) 



(41) 



where I is the short hand notation for /i, Pl,r = | (1 -F75) and q is the momentum transfer. 
fi's are the Wilson coefficients evaluated at the b quark mass scale in NLL order with values 



32| 



eff 



-0.31 , Cq = 4.154 , C 



10 



-4.261 . 



(42) 



The coefficient Cg" has a perturbative part and a resonance part which comes from the 
long distance effects due to the conversion of the real cc into the lepton pair . Hence, 
Cg" can be written as 

= Cg + Yis) + Cr , (43) 

where the function Y{s) denotes the perturbative par t coming from one loop matrix elements 
of the four quark operators and is given in Ref. 23|. The long distance resonance effect is 



given as 



33| 



cr = 4(3Ci + C2 + 3C3 + C4 + 3C5 + Ce) y ^ 



j/il),ip' 



rriy^ - s 



(44) 



where the phenomenological parameter n is taken to be 2.3, so as to reproduce the correct 
branching ratio B{B J/ipK* K*l+l-) = B{B J /i>K*)B{J /i) 1+1-). In this 
analysis, we will consider only the contributions arising from two dominant resonances i.e., 
J/ip and ip'. 

The matrix element for the decay Bg — )■ fi~^fi~'y can be obtained from that of the Bg — >■ 
H~^IJ,~ one by attaching the photon line to any of the charged external fermion lines. In order 
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to calculate the amplitude, when the photon is radiated from the initial fermions (structure 
dependent (SD) part), we need to evaluate the matrix elements of the quark currents present 
in (141 p between the emitted photon and the initial meson. These matrix elements can be 
obtained by considering the transition of a Bg meson to a virtual photon with momentum 
k. In this case the form factors depend on two variables, i.e., k"^ (the photon virtuality) and 
the square of momentum transfer = {pb — k)"^- By imposing gauge invariance, one can 
obtain several relations among the form factors at /c^ = 0. These relations can be used to 
reduce the number of independent form factors for the transition of the Bg meson to a real 
photon. Thus, the matrix elements for fi^ — t- 7 transition, induced by vector, axial-vector, 
tensor and pseudo-tensor currents can be parameterized as [34 1 



(7(fc,e)|s7^756|5,(pB)) = ie [e*^{pB ■ k) - {e* ■ pB)k^] 



Fa 



{^{k,e)\s^^h\Bs{pB)) 
{j{k,e)\sa^^q''-f5b\Bg{pB)) 
{^{k,e)\sa^,q''b\Bg{pB)) 



(^^fii/afS^ Pb k^ 



e [£*^iPB ■ k) - {e* ■ PB)kf,] Ft a 
ee^„ai3£*''p% k^Frv , 



(45) 



where e and k are the polarization vector and the four- momentum of photon, Pb is the 
momentum of initial Bg meson and Fj's are the various form factors. 
Thus, the matrix element describing the SD part takes the form 



M 



SD 



+ ^{e;{k ■ pb) - {e* ■ pB)k,,^ (^i 1^1 + B2 l^l^l 



(46) 



where 



Ai — 2Ct—t-Ftv + Cq- 



Ao = C 



—2C-j—Fta — Cg 



Fa 



10" 



mBs 



B, 



-C 



rriBs 



10" 



Fa 

^Bs 



(47) 



The form factors Fy and Fa have been calculated within the dispersion approach 35| . The 



dependence of the form factors are given as 



F{E,) = 13 



A + E^ 



(48) 
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TABLE II: The parameters for Bg ^ ^ form factors. 



Parameter 


Fv 




Fa 


Fta 


/3(GeV-i) 


0.28 


0.30 


0.26 


0.33 


A(GeV) 


0.04 


0.04 


0.30 


0.30 



where is the photon energy, which is related to the momentum transfer q"^ as 



1 - 



m 



Bs 



(49) 



The values of the parameters /3 and A are given in Table-2. The same ansatz (HHI) has 
also been assumed for the form factors Fta and Ftv- We use the decay constant of the Bg 
meson, which is evaluated in lattice QCD calculation as /^^ = 232 ± 10 MeV 36 1. 

When the photon is radiated from the outgoing lepton pairs, the internal bremsstrahlung 
(IB) part, the matrix element is given as 



M 



a 



IB 



VtbVtl fs, Cio 



p+ ■ k ■ k 



75 ^ 



(50) 



where p+ and p- are the momenta of emitted fi~^ and /i respectively. Thus, the total matrix 
element for the — ?> l^l^^ process is given as 



M = MsD + MiB ■ 



(51) 



The differential decay width of the Bg — ?■ fi^fi 7 process, in the rest frame of Bg meson is 
given as 



where 



(52) 



Ai = ^m|^(l-S)2i;,((S + 2n)(|Ai|2+|i?i|2) + (s-4n)(|A2p + |i?2 



64 



mo i - 



^C^?o((4n-^^-l)ln^ 



l + vi 



Vl 



2s Vl 



32 ri{l-rsf /bMCioAI 



(53) 



with s = q'^, s = s/m\^, ri = m^/m^^, vi = Jl — Am?^/q^. The physical region of s is 



4m^ < s < ml 
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The forward backward asymmetry is given as 



A 



FB 



Ai 



2m\;s{l - sfvf Re(^AlB2 + BIA2 



+ 32 /b, n(l-s)2ln(^)Re(Cio52 



(54) 



We have shown the variation of the differential decay distribution fl52p (in units of 10 
and the forward backward asymmetry (!54l) for Bs — > /i^/i~7 in Figure-4. 



fir- 
st 

.| 

■I 
:i 

'-\ 
■I 
1 -\ 



20 



0.5 



/ 
/ 

0.0 - — 



1 I 
-0.5 - \ 



/ \ 



25 30 5 10 15 



20 25 30 



FIG. 4: Variation of the differential branching ratio (in units of 10"^'') (left panel) and the forward- 
backward asymmetry with respect to the momentum transfer s (right panel) for the Bs — s- fi'^fj.^j 
process. 

As discussed earlier in the presence of fourth generation, the Wilson coefficients Cy g^io 
will be modified due to the new contributions arising from the virtual t' quark in the loop. 
Thus, these coefficients will be modified as 

Cr*(/i) = Cr(/i) + ^C^(/i), 



(55) 



The new coefficients C'-j 9 can be calculated at the My/ scale by replacing the t-quark mass 
by m[ in the loop functions as discussed in 23|]. These coefficients then to be evolved to 
the h scale using the the renormalization group equation. The values of the new Wilson 
coefficients at the m;, scale for m^/ = 400 GeV is given by C'-jirrih) = —0.375, Cg^rrib) = 5.831 
and C[q = -17.358. 
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Thus, one can obtain the differential branching ratio and the forward backward asym- 
metry in SM4 by replacing Cj^g^io in Eqs fl52|) and fl5^ by Cj°g iQ. Using the values of the 
Xf/ and (ps for m^/ = 400 GeV as discussed earlier, the differential branching ratio and the 
forward backward asymmetry for — t- fi^fi^'j is presented in Figure-5, where we have not 
considered the contributions from intermediate charmonium resonances. From the figure it 
can be seen that the differential branching ratio of this mode is significantly enhanced from 
its corresponding SM value whereas the forward backward asymmetry is slightly reduced 
with respect to its SM value. However, the zero-position of the FB asymmetry remains 
unchanged the fourth quark generation model. 

To obtain the branching ratios it is necessary to eliminate the background due to the 
resonances J/tpltp') with J/iplip') — > fj^^fj^'. We use the following veto windows to eliminate 
these backgrounds 

mj/^ - 0.02 < m^+^- < mj/^ + 0.02; 
m^i — 0.02 < m^+^- < m^i + 0.02. 

Furthermore, it should be noted that the |A^7bP has infrared singularity due to the emission 
of soft photon. Therefore, to obtain the branching ratio, we impose a cut on the photon 
energy, which will correspond to the experimental cut imposed on the minimum energy 
for the detectable photon. Requiring the photon energy to be larger than 25 MeV, i.e., 
> 5 uibJI, which corresponds to s < m\^{l — 5), and therefore, we set the cut 5 > 0.01. 
Thus, with the above defined veto windows and the infrared cutoff parameter, the total 
branching ratio for Bg — )■ fi^fi^'j process is found to be 

Br{Bs -> /^+/^"7) = 5.6 x 10"^ , (SM) 

< 2.2 X 10^^ , (for mt' = 400 GeV) . (56) 

The above branching ratio is comparable with that of the corresponding pure-leptonic pro- 
cess, Bs —7- /i^/U^, whose predicted branching ratio ^] for m^/ = 400 GeV is 

Br{B, fi+ii-) < 0.8 X 10"^ . (57) 



The LHCb 371] has searched for this process and set the upper limit as Bt{Bs — /i'''/i~) < 
1.2 (1.5) X 10"^ at 90% (95%) CL. Therefore, the Bs Ai+/i"7 decay channel could also be 
accessible there and hopefully it will be observed soon. 
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5 10 15 20 25 30 5 10 15 20 25 30 

s in GeV^ s in GeV^ 



FIG. 5: Variation of the differential branching ratio (left panel) and the forward-backward asym- 
metry with respect to the momentum transfer s (right panel) for the Bs — s- fi^fi^j process, in 
fourth quark generation model (red regions) whereas the corresponding SM values are shown in 
blue regions. 

VI. CONCLUSION 

In this paper we have studied some rare decays of the Bg meson in the fourth quark 
generation model. The large production of Bs mesons at the LHC opens up the possibility 
to study Bs meson with high statistical precision. The decay modes considered here are 
Bs —7- (f)n, Bs —7- 07, Bs — )■ 77 and Bs — )■ fi'^fi^-y, which are highly suppressed in the SM 
as they occurred only through one-loop diagrams. Therefore, they provide an ideal testing 
ground to look for new physics. The fourth generation model is a very simple extension of 
the SM with three generations and it can easily accommodate the observed anomalies in 
the B and Bs CP violation parameters for m^/ in the range of (400-600) GeV. We found 
that in the fourth generation model the branching ratios for these processes enhanced from 
their corresponding SM values. However, the mixing-induced CP asymmetry of in Bs — )■ (pir 
process enhanced significantly from its SM value. The CP violating observables in Bs — )■ 07 
are found to be small but nonzero. Some of these branching ratios are within the reach of 
LHCb experiments, hence the observation of these modes will provide us an indirect evidence 
for the existence of fourth quark generation. 
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